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ABSTRACT 

We measure the red giant branch bump (RGBB) of the Galactic bulge, the 
most metal-rich RGBB ever detected. The RGBB luminosity function peaks at 
the expected brightness, but its number density is very low relative to Galactic 
globular cluster calibrations, implying the Galactic bulge has a higher helium 
enrichment parameter AF/AZ > 4.0 for Y ~ 0.35 rather than the standard 2.0 
with Y = 0.27, which we suggest may be a common feature of galactic spheroids. 
The RGBB is (0.71 ± 0.02) mag fainter than the red clump (RC) in I toward 
the densest stellar regions imaged by the OGLE-III Galactic bulge photometric 
survey, (|/| < 4,2 < < 4). The number density of RGBB stars is (12.7 ± 2.0) 
% that of RC stars. The brightness dispersion of the RGBB is significantly lower 
than that of the RC, a result that is difficult to explain as the RGBB luminosity 
is known to significantly vary with metallicity. Sightlines that have two RCs have 
two RGBBs with similar properties to one another, an expected outcome if the 
Milky Way's bulge is X-shaped. We also find preliminary evidence of the Galactic 
bulge asymptotic giant branch bump, at a brightness of ~1.06 mag brighter than 
the RC in / and with a number density ~1.5% that of the RC. Accounting for 
the RGBB has a small effect on the best-fit parameters of the RC, shifting its 
best-fit peak brightness and reducing its brightness dispersion by ~0.015 mag 
each. 
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Introduction 



The red giant branch bump (RGBB) and the red clump (RC) are two observables of 
the color-magnitude diagrams (CMDs) of low and intermediate mass metal-rich populations 
which correspond to specific post main-sequence phases of stellar evolution. The RGBB 
occurs first, during ascent of the red giant (RG) branch. As the hydrogen burning shell 
expands, it eventually comes into contact with the convective envelope, which provides ad- 
ditional fuel. The star then gets fainter until the extra fuel is consumed whence it begins 
brightening again. As the star crosses the same luminosit y region three times, t here is an 
increase in the number density of stars at that luminosity ( ICassisi fc Salarislll997l ). The RC 
occurs after the degenerate helium core has accreted enough mass for core helium burning 
i nduced by the helium flash, it is the horizontal branch of an old, metal-rich population 
( IGirardi fc SalarisI l200l[ ). The relative number counts and brightness for these two popu- 
lations are a function of their com position. In particular, the RGBB gets fa i nter relative 



to the RC as metallicity i ncreases ( ICassisi &: SalarisI llQQTt IZoccali et al.l llQQQt iRiello et al. 



2003 : Di Cecco et al. 201o[). and the lifetime of the RGBB decreases with increasin g helium 
abundance (IBono et al.ll200ll : iDi Cecco et al.ll2010l ). while that of the RC increases (IRenzini 
1993). 



Whereas the RC has had a vibrant history as a trace r of Galact ic bulge reddening 
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2010f). there has been little study of the Galactic bulge RGBB. Indeed, we find only three 



brief mentions. IZoccali et al.l (|2003[ ) state that the RGBB is expected to be 0.7 mag fainter 
than the RC in J, and that the two will be vertically mixed on a CMD. The RGBB was also 
discussed in two recent investigations of the double RC toward sightlin e s near the bulge minor 
axis and at least 5 deg removed from the plane ( iNataf et al.l I2OIOI : iMc William fc Zoccali 
2OIOI ). Both works commented on a small, secondary excess of stars fainter than the RC, 
which they suggested may be due to the RGBB without further analysis. 



Interestingly, iMc William fc Zoccalil ( 120101 ) argued that the dou ble RC may have been 
observ ed before but not recognized as such. In their Figures 4 and 5, iBabusiaux &: Gilmore 
(l2005h plot a {J - K, K) CMD of a field at (/, h) = (0, 1) for which there is a second 
ov erdensity ^0.7 m a g fain ter than the RC. This second peak also shows up in an analysis 
by iNishiyama et al.l (120051 ) of photometry along the strip (|/| < 10,6 = +1). We argue the 



secondary density was in fact the RGBB and not the double RC. Within our own data, we 
find that the RGBB signal is stronger for sightlines closer to the plane, an effect that can be 
explained by the lower geometric (and thus brightness) dispersion for those sightlines. 



5x104 




-1 







1 



Fig. 1. — The magnitude distribution for 2,836,575 OGLE-III point sources that are within 
4 deg^ of (/, b) = (0, 0), that have an / measurement within 1.5 mag of the local RC centroid, 
and that are no more than 0.35 mag bluer than the nearest RC. The top panel shows the fit 
to a Gaussian plus a power-law for the RC and RC sources, the middle panel shows the fit 
once a second Gaussian for the RGBB is added. Ax^ = 6670. The bottom panel shows the 
residual to the latter fit. Prominent and significant are a few non-Gaussianities to the RC 
and an excess of stars ~1.1 mag brighter than the RC, hkely due to the AGBB. 
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In this Letter, we use the OGLE-III Galactic bulge VI photometry dataset and find 
that the RGBB is ~0.71 mag fainter than the RC in J, and that its number density is 
(12.7 ± 2.0)% that of the RC. Accounting for the RGBB slightly affects the parameters of 
the RC, the number density of RC stars shrinks by ~2%, the peak brightness shifts by ~0.015 
mag toward brighter mag, and the brightness dispersion decreases by ~0.015 mag. The fit 
and distribution is shown in Figure [H We may also have detected the asymptotic giant 
branch bump (AGBB), at ~1.1 mag brighter than the RC with a number density of order 
~1.5% that of the RC. Sightlines with two RCs have two RGBBs with similar properties. 
In section [2|, we describe the OGLE-III observations. We summarize the evidence for the 
Galactic bulge RGBB in section [31 We present the method of the clump-centric CMD in 
section HI and our method of measurement in section [51 Each of the measured RGBB 
parameters are presented and interpreted in section [SI We discuss our findings as well as the 
prospects and directions of future research in section [HI 
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Fig. 2. — The OGLE-III Galactic bulge observing coverage, color-coded by the {V — I) color 
of the RC. The RC gets redder closer to the plane due to higher interstellar e xtinction and 
redde ning - the intrinsic color of the Galactic bulge RC is {V — 1) ~1.08 ( iBensby et al. 
20101). As the RGBB is ~0.71 mag fainter than the RC, and OGLE-III is complete down to 



21st magnitude, we are able to detect the RGBB over nearly the entire field of view. 
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Data 



OGLE-III observations were obtained from the 1.3 meter Warsaw Telescope, located at 
the Las Campanas Observatory in Chile. The camera has eight 2048x4096 detectors, with a 
combined field of view of 0.6° x 0.6° yielding a scale of approximately 0.26"/pixel. We use 
observations from a mosaic of 267 fields directed toward the Galactic Bulge, which are almost 
entirely within the range < |/| < 10 and 2 < |6| < 7. The photometric coverage is shown in 
Figure |2j Of the 267 fields, 37 are at northern latitudes. In this study we do not investigate 
the fields BLG200, 201, 202 and 203, located in and around {l,b) = (-11,-2.5), due to 
the anomalous CMDs toward those sightlines, possibly due to a smaller scale differential 
reddening. More detailed descriptions of the in s trume ntati on, photometr i c redu ctions and 
astrometric calibrations are available in lUdalskil ( 120031 ) and lUdalski et al.l (|2008[ ). 



Our quantitative results are obtained from 2,836,575 OGLE-III point sources that are 
in subfields centered within 4 deg^ of (/, 6) = (0,0). We chose this subsample as it is 
more robust - it has a higher density of stars, lower disk contamination, lower geometric 
dispersion, no double clump, and no controversy of multiple Galactic bars implying distinct 
stellar populations. Our error estimates are based on comparisons with other coordinate 
selection criteria. 



3. Three Lines of Evidence Toward the Red Giant Branch Bump 

There are three lines of evidence that lead us to the conclusion that the RGBB is a 
measurable feature of Galactic bulge (/, V — I) CMDs. It appears as a statistically signifi- 
cant secondary peak in our magnitude histograms, there are two such additional peaks for 
sightlines with two RCs, and the peak is approximately where it is expected to be from 
globular cluster calibrations. 

The first can be recognized from the distribution of differences between the observed 
distribution of RC-I-RG stars and their model fit of a power-law plus a Gaussian. There is a 
clear excess above the fit for stars 0.5-1.0 mag fainter than the RC. The absolute area of the 
residuals is reduced by up to 75% once a second Gaussian is added to account for the RGBB. 
XdoF ^he fit is 75.38 without an RGBB. When a free Gaussian for the RGBB are added, 
XdoF drops to 5.34 - an unambiguous statistical signal. We note that the large remaining 
value of XdoF indicates that 8 free parameters are not sufficient to fully describe the data. 
That is not surprising as a variable slope to the red giant branch, differential reddening on 
small scales, foreground disk contamination, an asymptotic giant branch bump, and non- 
Gaussianities in the RC and RGBB are all legitimate possibilities for additional parameters. 



- 6 - 



0.5 |— • Galactic Globular Clusters 
Dwarf Spheroidals 



m 

om 

> 
< 



-0.5 - 



-1 



T 1 1 1 1 1 1 1 1 1 1 1 1 r 

• Milky Way Bulge 
□ cj Cen 




I iV 




-2 



-1 
[M/H] 







Fig. 3. — The ^Vzahb fo'^ the 15 Galactic globular clusters analyzed by iDi Cecco et al. 



(120 lOl ) are shown as red circles, the 7 RGBBs of the cluster uj Cen iRey et al.l (120041 ) are shown 
as empty blue squares, and the 12 RGBB measurements toward 9 local dwarf spheroidals are 
shown as filled magenta triangles (Monachesi et al. 2011, Monelli et al. 2010, and references 
therein) . 



The second line of evidence is related to the double-RC of the Galactic bulge, seen on 
sightlines within ~2 degrees of the bulge minor axis and at least ~5 degrees removed from 
the plane. The two RCs are approximately equally populated, have equal or very nearly 
equal {V — I) and (J — K) color, and are separated in brightness by ~0.5 mag ( iNataf et al. 



2010l : iMcWilliam fc Zoccalill2010l ). For double-clump fields, combining data from different 
sightlines with different reddening is non-trivial, however we can use the residuals to a single- 
clump fit on small scales as a preliminary diagnostic. The signal is clear: two clumps and two 
bumps. If the two clump-bump pairs are found to have the same separation in brightness, 
and the same population fraction, it would rule out explanations of the double clump based 
on age and chemistry and confirm an X-shaped Milky Way bulge as the explanation. Both 
those statements appear approximately correct based on the distribution of residuals. If 
detailed analysis confirm that the two populations have the same metallicity as one another. 
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they will also show whether they have the same metallicity as the bulge giants toward 
sightlines closer to the plane, th ereby providiii g an in dependent test and calibration of the 
metallicity gradient detected by IZoccali et al.l (|2008[ ). Unfortunately, the method used in 
most of this paper cannot be directly applied to double-clump fields, due to the varying 
separation between the two RCs and their varying relative population densities, as well as 
the fact these variations are not yet well-parametrized. In Figure H] we plot the sum of 
residuals to the fits of single RC + KG distribution to thousands of sightlines toward both 
single-clump fields and double-clump fields. The residuals for the double-clump fields betray 
not only the presence of a second RC, but also of two RGBBs. 



The third line of evidence is that the brightness peak of the RGBB we report here, 
0.71, is consistent with that expected from stellar evolution models and globular 
cluster calib r ation s. In their analysis of 54 Galactic globular clusters observed with HST, 
Riello et al.l ( l2003l ) reported a difference in brightness AF555W between the RGBB and 
the Zero Age Horizontal Branch (ZAHB) evaluated at the RR Lyrae instability strip. They 
found that A F555W = 0.45 for NGC5927 and 0.53 for NGC6624, two metal-rich globu- 



lar clusters (Carretta et al. 



20091). For the me t alhcity of the Gala ctic bulge, [M/H] ^ 0.0 



(iFulbright et al.l l2006l. 120071: Zoccah et al.ll2008l : iBensbv et al.ll2010f l. the best-fit model pre- 
sented by iRiello et al.l ( l2003l ) has AF555H^ = 0.639 mag for a 10 Gyr stellar population. 



and A.F555Vt^ = 0.74 1 mag for a 12 Gyr stellar population. Similar results are obtained by 
Di Cecco et al.l ( l2010l ) in an investigation of AVzahb 15 Galactic globular clusters. We 



plot various examples from the literature as well as our own in Figure Ei 



4. A New Diagnostic for Stellar Populations: The Clump-Centric 

Color-Magnitude Diagram 

We introduce the clump-centric color-magnitude diagram (CCCMD) as a diagnostic 
means to study the fine details of large stellar populations in the presence of differential 
reddening and geometry. It is the color and magnitude of every star within a large angular 
range relative to the nearest RC centroid. 

CMDs are among the most powerful observational tools in Astronomy. Unfortunately, 
optimizing their construction for stellar populations such as that of the Galactic bulge is a 
non-trivial task. Over small angular scales, their diagnostic power will be limited by Poisson 
noise, restricting the analysis of a stellar population to only its most promin ent features 



For e xample, within the CMDs used in our previous work on the double RC ( iNataf et al. 



2010|), we used CMDs just large enough to contain 1200 RC-I-RG stars, so as to minimize any 



possible gradients. With that number we would only expect ~30 RGBB stars per sightline. 
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rendering it a difficult population to identify let alone classify as it is mixed with the RC 
and RG populations. 
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Fig. 4. — The sum of a residuals to a single-chimp+power-Iaw fit over thousands of sightlines 
containing no more than ~1200 RC+RG stars to two general areas the bulge, (|/|<1, |&|<3) 
in the top panel, and (|/|<1,— 8<6< — 5.5) in the bottom panel. The residuals in the 
top panel, toward single-clump areas, are dominated by the RGBB and the requirement of 
normahzation. Those toward the bottom panel reveal the presence of two RGs, as well as 
two RGBBs that are approximately equally well-populated and equally well-separated to 
their respective RGs. 
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O yer la rge angular scales, differential reddening and geometry both become problems. 

Sumil ( 120041 ) found that the optical reddening E{V — I) can vary by over ~1 mag within 
1 degree. Worse, the extent of differential reddening is larger for sightlines closer to the 
plane, which have a higher density of stars. This would produce a diagonal smearing over 
any large CMD. This smearing would have no clear form as both the slope of the reddening 
law, Rj = All E{y — I) and the total reddening have patchy and discontinuous variations. 
Dere ddening is typical l y used as a solution to this problem. In their study of the Galactic 
bar, iRattenbury et al.l ( l2007al ) dereddened each field to the intrinsic color of the RC. An 
issue with this method is that the slope of the reddening law is known to vary with direction 
even though it is not known precisely how it varies with direction. A change in Ri of ~0.1 
over 2 magnitudes of reddening leads to an error of 0.2 mag in the dereddened brightness. 

Further, it is also known that within the coordinate range (—4 < / < 4) the brightness 
of G alactic bulge RC stars varies by ~ 0.4 mag due to the orientation of the Galactic 



bar (jStanek et al.l 119971 : iBabusiaux fc Gilmord l2005l : IRattenbury et al.l |2007aJ) . This effect 
is in some ways more fundamental than the first effect, as it cannot be evaded by using 
"reddening-free" near-IR or mid-IR data. The obvious countermeasure to assume a Galactic 
bar model would be an imperfect solution. The actual offset due to the bar is unknown at 
the level of ~0.1 mag. In a ddition, the bright r iess v ariation with longitude is known not to 
be independent of latitude. iNishiyama et al.l ( 120051 ) found that the variation of brightness 
with longitude for sightlines toward 6 = 1 is shallower than that found further from the 
plane, an effect that could be due to a additional inner bar. 

These two factors can be almost completely eliminated by using a CCCMD. By con- 
structing a grid of ~30,000 RC {V — J, I) measurements across the OGLE-III bulge sky, we 
achieve no worse than 4.5' resolution. Over that range, the variation in geometry is negli- 
gible, and the effect of differential reddening is minimized to its expected extent over that 
length scale - typically a few hundredths of a magnitude or less. By amassing a larger number 
of stars, we can make detailed measurements of structures such as the RGBB and possibly 
the asymptotic giant branch bump (AGBB), as well as investigate possible systematics in 
the RC parameters. 



5. Fitting for the Red Giant Branch Bump 



We fit for the properties of the RGBB in a two-step pr ocess. We fi r st fit for the RC 
without accounting for the RGBB in a manner described in iNataf et al.l ( 120101 ). The CC- 
CMD is produced by subtracting from each star the {V — I, I) color and magnitude of the 
nearest measured RC. Stars that are nearest to the RCs located in zones of rapid differential 
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reddening, low stellar density, or Vrc > 19.5 are left out of the CCCMD. 



We then fit for the combined RC+RG+RGBB population using the clump-centric color- 
magnitude cut 

- 0.35 <{V- I)cc 
-1.5 < Ice < 1.5, (1) 



and the model 



N{'m) = A exp 



B(I-I 



RC) 



RC 



"^rrciRc 



exp 



RC) 



24c 



N, 



RGBB 



'^T^CFRGBB 



exp 



RGBB, 
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where A, B, Nrc, Irc, and aRc are as described in iNataf et al.l (|2010[ ): and Nrcbb, Irgbb, 
and ctrcbb are the analogous parameters for the RGBB. 



6. Four Measurable Parameters for the RGBB 

The fitting routine of section \5\ leads to four measurable parameters for the RGBB: 

Irgbb — Irc (3) 

(4) 

(5) 
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6.1. The Brightness Parameter AI^^^^ 

The first parameter, AI^'^^^ , is the difference in peak brightness between the RGBB 
and RC, and is 0.71±0.02 mag. This is consistent with the expectation from globular clusters. 
We use / rather than V because the RC has a thinner Gaussian in J, allowing a clearer 
separa tion of the two populations. This is expected. In the models of iGirardi fc Safaris 
teOOll ). My^Rc for an old population has variati ons with metallicity ~2x that of Mj^rc, an 
effect confirmed in the empirical investigation of iPietrzyhski et al.l ( l2010l ). who looked at the 
RC of 15 nearby galaxies observed with HST. Moreover, any residual differential reddening 
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will be ~2x as significant in V. These two effects render the bulge RC non-horizontal in V, 
further complicating the fitting routine. 

There is a difference between this parameter and that predominantly used in the liter- 
ature. We computed AI^^^^ whereas most results present /SV^ahb ■ However, these two 
values should be very nearly equal as the two largest biases are not large and go in opposing 
directions. Firstly, AV^^^^ should be a little larger than AI^q^^ because the RC will be 
a little bluer than the RG stars, the bias is expected from stellar theory but is consistent 
with negligible in our data. Conversely, AV^^^^ will be a little smaller than AV^^^^ as 
the ZAHB is the dimmest phase of horizontal branch evolution, however in our analysis of 
globular cluster data we find that this effect can be no more than ~0.05 mag. We thus adopt 
the approximation AI§^^^ = AV^^§^ . 

It is reassuring that the Galactic bulge has the faintest RGBB relative to its horizontal 
branch as it is the most metal-rich RGBB detected thus far. Galactic globular clusters do 
not typically reach met allicities as high as [ M/H] ^ 0.0, and those that do have substantial 



differ ential reddening (lOrtolani et al.l 120011 ) or multiple stellar populations (iFerraro et al. 



20091 ). effects that render the RGBB harder to detect. This expansion of the parameter space 
at the metal rich end follows recent, complementary detections of the RGBB in younger and 
more metal-poor systems, those of the nearby dwarf gala xies. There have beeii detections of 
the RGBB toward the Sculptor dwa rf spheroidal galaxy (IMajewski et al.lll999l). Ursa M inor 
(IBellazzini et al.ll2002l). Sagi tt arius (IMonaco et al.ll2002[) . and Sextans (iLee et al.ll2003l ). and 
M32 (IMonachesi et al.ll201ll ). iMonelli et al.l (120101 ) recently reported on the detection of the 



RGBB toward Cetus, IC1613, LGS 3 and Tucana. The observed relation between AVzahb 



and [M/H] is shown in Figure [3], with the IZinn fc WestI (119841 ) metallicity scale assumed for 
the Galactic globular clusters and u Cen. 



We estimate [M/H] via the conversion function suggested by ISalaris et al.l (119931 ): 

[M/H] = [Fe/H] + log(0.638 * IQi^/^^l + 0.362) (7) 

and an [a/Fe]= -1-0.4 for the Galactic globular clusters and u Cen, and [a/Fe]= +0.25 for 
the Galactic bulge and M32. 



6.2. The Population Parameters fsump and EWrgbb 

The total number density of Galactic bulge RGBB stars is down by a factor of ~2-3 
from that expected from globular cluster calibrations, either relative to the RC or to the RG 
branch. For both parameters, the Galactic bulge stellar population is a lot more similar to 
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that of M32 (iMonachesi et al.ll201ll ) than to that of Galactic globular clusters at comparable 
metallicity, which suggests a common factor in the stellar evolution of galactic spheroids. 



We find that the ratio of the number of RGBB stars to RC stars is / = (12.7 ± 2.uj7o, 
and that the equivalent width of the RGBB EWrgbb = (0.104 ± 0.020) mag. Both of them 
parametrize the total population of the RGBB, but the former does so by comparing its num- 
bers to that of the RC, and the latter does so relative to that of the underlying RG branch. 
Both have their relative strengths, / is less likely to be biased by disk-contamination, whereas 
EWrgbb may be easier to interpret theoretically as it compares two p opulations before the 
helium flash. We do not use the standard strength calibrator RBump (IBono et al.ll200ll ) as 
it would be extremely vulnerable to disk contamination and photometric incompleteness for 
stars fainter than ^ ^ 21. 

As these are both new parameters previously undefined in the literature, we had to 
compute their value in globular clusters. We first used a s ubsample of globular clusters 
observ ed with the Hubble S pace Telescope (jPiotto et al.ll2002l ). that have an updated metal- 
licity (ICarretta et al.l 120091 ) . and are densely-populated. We correct for their photometric 
completeness functions. This gave us the value of EWrgbb for 30 globular clusters. The 
mean and sample deviation are 0.32 and 0.18 mag, and there is a statistically significant 
correlation with metalhcity, EWrgbb = 0.370 + 0.133[Fe/H]. We then took the subsample 
of 7 metal-rich globular clusters (NGC104, 5927, 6356, 6624, 6569, 6637, and 6760) that are 
well-sampled, do not suffer from differential reddening and for which all or nearly all hori- 
zontal branch stars are RHBs rather than BHBs. The mean and sample deviation of their 
fractions are 31.7% and 9.7%. The correlation with metallicity is / = (22.2 — 10.7[Fe/H])%. 
However there are only 7 points spanning a small range of metallicities relative to their 
extrapolation. Both sets of results are shown in Table [H 

The result of low number density should not be surprising. If indeed the RGBB should 
be more easily measured in metal-rich populations, why was that of the Galactic bulge not 
measured decades ago? We argue that if the bulge RGBB had 20-30% the number density 
of the RC it would have likely been detected, interpreted and integrated into the literature 
by this time. 

We suggest that the Galactic bulge has a higher helium to metallicity enrichment param- 
eter (AY/AZ) than Galactic globular clusters, typically assumed to be ~2. It is expected 
from stellar evolu tion models that helium enrichment reduces the lifetime of the RGBB. 
Bono et al.l (120011 ) showed that for a 14 Gyr population, Rbump should decrea se by ~10%, for 



an inc rease in the initial stellar helium abundance AY = 0.06. More recently. iDi Cecco et al. 



(120 lOl ) estimated that for an [M/H]= —0.35 population, the lifetime of the RGBB will be 
decreased by one half as Y is increased from 0.26 to 0.31, and by two thirds when Y = 0.33. 
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The numbers will be a little higher for a [M/H] ^ population. An increase d abundance o f 
helium is also expected to lead to a longer lifetime for the horizontal branch (jRenzinilll994j ). 
making fBump an even further sensitive parameter to helium. 
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Fig. 5. — TOP: The RGBB:RC ratio for 7 metal-rich globular clusters are shown as triangles, 
and that of the Galactic bulge is shown as the pentagon. BOTTOM: EWrgbb for 30 globular 
clusters are shown as triangles, and that of the Galactic bulge is shown as the pentagon. The 
solid black line is the theoretical relationship EWrgbb = 0.404 + 0.134[Fe/ H] for a 10 Gyr 



popul ation. For both panels, the globular c luster metallicit y is ta ken from iCarretta et al. 
(120091 ) ■ and that of the Galactic bulge from iFulbright et al.l ( 120061 ). The data for the local 
dwarf spheroidal M32 is shown as the empty black square and is taken from lMonachesi et al. 
(j2011 ). for which we assume a metallicity error of 0.1 decs. 



We have estimat ed the theoretical relationships using the YREC evolutionary code 
(jPelahaye et al.ll2010l ). We assumed a primordial helium abundance of 25% and a solar 
helium abundance of 27% with a linear enrichment parameter AF/AZ in between, and we 
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find EWrgbb = 0.404 + 0.134[Fe/H]. For a [Fe/H]= star, the theoretical lifetime of the 
RGBB phase is tjicBB = 23.6 — 1.4(100F — 27) Myr, as suming a 10 G yr population for both 
calculations. Since the lifetime of the RC is ~100 Myr (lYi et al.lll997l ). this implies an initial 
helium content of ~35% by mass for the typical metallicities of bulge stars. A higher helium 
abundance would be required to fit for the result EWrgbb = (0.104 ± 0.020) mag, but that 
parameter is significantly more sensitive to age effects and disk contamination. We note that 
EWrgbb decreases with age, whereas ftump increases with age. The theoretical relation at 
[Fe/H]= and t= 10 Gyr is EWrgbb ~ 0.408 - 0.024(100F - 27) - 0.008(t - lOGyr) mag, 
requiring an initial helium content of 39%. We report 35% as our initial estimate, but we 
recognize that a more concrete estimate must account for the impact of disk contamination 
as well as integrating over a well-motivated composite population in metallicity and helium 
enrichment, tasks beyond the scope of this work. The effect of helium enrichment on the 
lifetime of the RGBB phase for a 1 Mq star with solar metallicity is shown in Figure O 



6.3. The Dispersion Parameter Aa^ 



The best-fit dispersion in brightness of the RGBB stars is lower than that of the RC 
stars. For our canonical field, we find ctrg = 0.253 ± 0.012, and ctrgbb = 0.203 ± 0.003, 
implying a difference in dispersion Ao"^ = 0.05, or Aa = 0.15. This is surprising as the 
luminosity of the RGBB is a steeper function of metallicity than that of the RC. 

We estimate the expected intrinsic width of the RGB B as follows. We fir st take 195 
RG+RC [Fe/H] measurements toward Baade's wind ow from IZoccali et al.l (j2008[ ). which we 
then convert to [M/H] using the same prescription as lBrown et al.l ( l2010l ): [a/Fe]= -1-0.30 for 
[Fe/H] < 0.0, [Qi/Fe] = 0.0 for [F e /H] > 0.5, and interpolat i ng liri early in between. We then 
combine the results of iRev et all koo4 ] and lPi Cecco et al.l koid ). getting AI^^^^ = 0A7+ 
0.53[M/H], assuming AI^q^^ = AV^'ahb ^ before and thereby allowing us to empirically 
estimate the brightness of each metallicity population. However, the brightness distribution 
is not just a linear transform of the metallicity distribution, as more metal-rich RGBB stars 
go through that phase earlier in their RG ascent, and thus the RGBB distribution will be 
biased toward the more metal-rich subset of stars. This effect is the same reason why the 
RGBB is harder to detect and classify in metal-poor clusters, in which the RGBB is brighter. 
We estimate the size of this bias by weighting each metallicity point by the observed relative 
number density of the Galactic bulge RG branch at that brightness, obtaining the following 
characteristics of RGBB stars: 



0.47 + 0.53[M//f] 



Wi 



exp 



(8) 
(9) 
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11.94 11.96 11.98 12 9.68 9.7 9.72 9.74 

Age (Gyr) Age (Gyr) 




Age (Gyr) Age (Gyr) 

Fig. 6. — The theoretical RGBB for 1 Mq stars with solar metallicity, but with an initial 
helium content of 27,30,33 and 36%, going clockwise from the top left. Blue dots represent 
the phase of RG evolution at the luminosity of the RGBB, and red dots the adjoining ascent 
of the red giant branch. The green lines are the tangent lines to the luminosity evolution 
before and after the RGBB phase, from which the lifetime is estimated. As the RG lifetime 
is very short compared to the total stellar lifetime, the luminosity evolution of a single star 
will be an excellent approximation to the luminosity function of a large number of stars with 
very small variations in initial mass. 
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where /j is the estimated clump-centric brightness of the RGBB for that metalhcity, 
Wi is the weight of that brightness component in the estimated distribution, and B is the 
measured exponential slope of the RG branch, 0.547. We calculate the weighted mean and 
variance: 

7=^J^ (10) 

= ^' > (11) 



obtaining a mean J = 0.560 and ct/ = 0.157. This is reassuringly lower than the measured 
dispersion ajiGBB = 0.203 mag, a "requirement" as there will be additional variation induced 
by differential reddening, binary blending contamination and distance variation due to the 
large physical size of the bulge - contributions expected to similarly affect the RC. We 
briefly note the difference in the mean compared to our measurement could also be removed 
by discarding the small sample of expected RGBB (less than 5%) that are within 0.25 mag 
of the RC as our fitting routine is not sensitive to those RGBB stars - the mode of our 
expected distribution is at ~0.71 mag. 

A smaller brightness dispersion of the RGBB than for the RC is a surprise. Performing 
the same calculation as above for RC stars, assuming li fetimes independent of m etallicity and 



a brightness relation A//j(7/A[Fe/H] = 0.14 mag/dec (jPietrzyhski et al.l 120101 ) implies a dis 



persion of o"/ ~ 0.06 mag. Galactic bulge RG sta rs have been measured to h ave similar proper 



motion distributions as their RC counterparts ( iRattenbury et al.ll2007bl l. an observational 
constraint that negates the possibility of a substantial bias- inducing age-metallicity-geometry 
correlation. There could be a bias induced in the parameters if the bulge RC has a skewed 
distribution as it is fit to a Gaussian, the residuals would then affect the fit to the RGBB. 
We do not expect this to be the case as the brightness distribution of the RC should be made 
very symmetric by geometrical dispersion. Forcing the intrinsic dispersion of the RGBB to 
be 0.05 mag higher than that of the RC increases by ~320, demonstrating that a signif- 
icantly different functional form would need to be assumed and justified to obtain a larger 
dispersion for the RGBB. We therefore conclude that the narrow luminosity function for the 
RGBB relative to the RC may be real. This could be achieved if there are subpopulations 
in the bulge that contribute substantially to the RC population and thus its dispersion but 
very weakly to the RGBB, such as a much more helium-rich population. In this scenario, it 
is the RC dispersion is too high rather than the RGBB dispersion being too low. 



-17- 



7. The Asymptotic Giant Branch Bump? 



The AGBB (lGallartlll998l ) may be present as a second peak in our magnitude distribu- 
tion, at a brightness ~1.06 mag higher than that of the RC, visible in the bottom panel of 
Figure [TJ Since the AGBB is a weak signal, it is necessary to constrain a third Gaussian in 
order to properly fit it, and we first do this by applying a penahy to solutions with a peak 
AGBB brightness outside the interval 1.00 < Irc — Irgbb < 1-30. The best-fit then has a 

~ 125, with a peak-brightness ~1.12 mag brighter than the RC and a population of 
0.8% that of the RC. However, it has a seemingly unphysical dispersion of ~0.11 mag, which 
is substantially narrower than that of the RC and the RGBB. This could be the case if the 
AGBB probes a much narrower range of stellar evolution than the RC or RGBB or if it is 
more sensitivity to an age-geometry-metallicity correlation, but we suspect it is like due to 
systematics in the fit. Constraining the AGBB to have the same dispersion as the RC yields 
a solution where the AGBB is ~1.04 mag brighter than the RC in /, has a number density 
that is 1.8% that of the RC, with a Ax^ ~ 79. Constraining it to have the same dispersion 
as the RGBB yields an AGBB to RC ratio of ~1.5%, a separation in brightness of ~1.08 
mag, and A^^ ~ 100. The signal for the AGBB is not due to the RC of a foreground spiral 
arm or Galactic ring as the separation in brightness is independent of longitude, ruling out 
different geometrical configurations for the two peaks. It also shows up in the double-clump 
fields as seen in Figure IH This feature could provide an important insight into Galactic 
bulge evolution in the future should theoretical investigations and empirical calibrations 
better constrain its detailed properties. We note that the Galactic bulge A GBB:RC ratio of 



1.5 ± 0.5)% is substantially lower than the (6.0 ± 1.8)% detected in M32 ( iMonachesi et al. 



201 ll ). which is perhaps surprising given that their RGBBs behave in a similar manner. 



8. Discussion and Conclusion 

We have demonstrated that the Galactic bulge RGBB peaks at or near the expected 
luminosity, but that its number-density is several times lower than would be expected from 
globular cluster calibrations. The brightness dispersion of the RGBB is lower than that of 
the RC. Both effects can be explained if there is a substantial, helium-rich population in the 
bulge that contributes to the RC but not the RGBB. 

This work will be difficult to extend to most of the Galactic bulge as many sightlines 
have much higher brightness dispersions due to geometrical effects, leading to degeneracies 
in the fit to the RC and the RGBB. However, there are several promising avenues for further 
investigation. With our dataset, investigations of the AGBB and of the double clump-bump 
pairs toward what are probably the X-wings of the Milky Way bulge are both possible. 
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Additionally, any large photometric dataset toward the inner few square degrees will be 
probing sightlines with minimal geometric dispersion, allowing a clearer separation of the 
RC and the RGBB. These are also the densest stellar fields, and so it might be possible to 
measure gradients. 



We find one literature mention that the Galactic bulge may be helium-rich. iRenzini 



(119941 ) estimated 2 < AF/AZ < 3, for a helium content 0.30 < Y < 0.35. His analysis was 
based on the high value of the number ratio Rc of R C to RG stars, an effect previously argued 
to be due to a relatively young stellar population ( jPaczynski et al.lll994( ). We comment on 
one chemical-enrichment model of interest. In their analys i s of th e split horizontal branch 
of the metal-rich globular cluster Terzan 5, iD'Antona et al.l ( 120101 ) showed that the helium- 
rich ejecta of massive asymptotic giant branch stars could lead to a substantially helium- 
enriched (Ay ~ 0.07) secondary population within ~ 100 Myr. T his is shorter than the 
currently expected ^0.5 Gyr timescale of Galactic bulge formation (ICescutti fc Matteucci 



2010l : iBrown et al.l 120101 ). As our method generalizes well to any large stellar system, we 
wonder if a similarly low RGBB population density is to be found in the pseudobulges of 
nearby massive spiral galaxies, and if enhanced helium-enrichment is a characteristic of such 
populations. 

Enhanced helium enrichment, if confirmed, will lead to a downward revision of the 10-f 



Gyr age of the Galactic bulge determined from photometric isochrone fitting ( iZoccali et al. 



20031 : iBrown et al.ll2010l ). This would possibly bring it in line \ yith the substantial ly younger 
ages detected via microlensing spectroscopy of dwarf stars ( iBensby et al.l l2010l ) and star 
counts of infrared carbon stars ( ICole fc Weinbergll2002l ). 
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Table 1: Measured globular cluster population parameters for the HS T sample of iPiotto et al 
( 2002 ). [Fe/H] is taken from the globular cluster metallicity scale of Carretta et al.l ( 2009 ). 



Object 


Nrgbb 


Nrc 


EWrgbb 


[Fe/H] 


[Fe/H] 


NGC 104 


87 


306 


0.35 


-0.76 


0.02 


NGC 362 


38 


- 


0.36 


-1.30 


0.04 


NGC 1261 


14 


- 


0.20 


-1.27 


0.08 


NGC 1851 


42 


- 


0.39 


-1.18 


0.08 


NGC 1904 


11 


- 


0.14 


-1.58 


0.02 


NGC 2808 


121 


- 


0.38 


-1.18 


0.04 


NGC 5024 


14 


- 


0.16 


-2.06 


0.09 


NGC 5824 


18 


- 


0.07 


-1.94 


0.14 


NGC 5904 


32 


- 


0.61 


-1.33 


0.02 


NGC 5927 


38 


176 


0.38 


-0.29 


0.07 


NGC 5986 


19 


- 


0.20 


-1.63 


0.08 


NGC 6093 


23 
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0.18 


-1.75 


0.08 


NGC 6139 


36 


- 


0.29 


-1.71 


0.09 


NGC 6229 


24 


- 


0.16 


-1.43 


0.09 


NGC 6266 


18 


- 


0.07 


-1.18 


0.07 


NGC 6284 


31 


- 


0.57 


-1.31 


0.09 


NGC 6356 


88 


364 


0.37 


-0.35 


0.14 


NGC 6388 


188 




0.28 


-0.45 


0.04 


NGC 6402 


51 




0.42 


-1.39 


0.09 


NGC 6441 


257 




0.31 


-0.44 


0.07 


NGC 6522 


46 




0.77 


-1.45 


0.08 


NGC 6569 


35 
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0.28 


-0.72 


0.14 


NGC 6624 


42 


83 


0.44 


-0.42 


0.07 


NGC 6637 


45 
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0.76 


-0.59 


0.07 


NGC 6638 
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0.18 


-0.99 


0.07 


NGC 6760 


47 


139 


0.37 


-0.40 


0.14 


NGC 6864 
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0.32 


-1.29 


0.14 


NGC 6934 


18 




0.40 


-1.56 


0.09 


NGC 7089 
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0.07 
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0.23 


-0.37 


0.07 



